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FEAP74,  a  finite  element  code  developed  for  USDOT  to  treat  contact-impact  problems,  was  secured  for 
the  analytical  portion.  Viscoelastic  and  large  deflection  elements  in  the  code  were  defective,  and  analysis 
was  limited  to  elastic  impact.  The  program  was  used  to  analyze  a  considerable  number  of  long,  laminated 
plates  clamped  on  two  parallel  edges.  The  analysis  treated  flexure,  contact,  and  transverse  shear  defor¬ 
mations  together  and  separately  to  obtain  the  relative  orders  of  magnitude  of  each  effect. 

Tests  were  designed  to  duplicate  the  conditions  of  the  analysis.  AS-3501  graphite/epoxy  plate  samples 
were  fabricated  and  tested  using  the  NAVAIRDEVCEN  impact  tower  with  automated  data  storage  and 
manipulation  capabilities.  Each  test  specimen  was  strain  gaged  so  that  time-histories  of  back-surface  strain 
could  be  measured  and  compared  with  analysis. 

The  program  has  shown  that  flexural  viscoelasticity  is  negligible,  but  that  viscoelastic  transverse  shear 
and  contact  deformations  are  important  for  short  composite  plates.  It  appears  that  membrane  effects  may  be 
important  for  longer  plates  where  strains  are  significantly  underpredicted  by  elastic  small-deflection  analysis. 
These  results  will  give  guidance  to  analysts  who  must  predict  strains  and  stresses  in  order  to  achieve  safe, 
efficient,  impact-resistant  composite  structures. 


Since  previous  creep  data  on  AS-3501  graphite/epoxy  material  was  determined  using  minutes  as  data 
intervals,  it  was  unsuitable  for  time  frames  of  a  fraction  of  a  second  which  represent  duration  of  impact  load¬ 
ing.  The  determination  under  this  program  of  an  experimental  viscoelastic  creep  compliance  for  AS-3501 
graphite/epoxy  material  will  provide  analytical  viscoelastic  models  with  accurate  property  data  with  which  to 
analyze  laminates  undergoing  low-velocity  impact. 
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I.  INTRODUCTION 


Laminated  fiber  composites,  used  extensively  in  military 
aircraft,  can  be  damaged  by  impact  of  blunt,  hard  objects  such  as 
stones  or  dropped  tools.  It  is  necessary  to  be  able  to  analyze 
these  structures  and  design  them  to  resist  impact  damage.  For  low 
velocity  impact  where  wave  propagation  may  not  play  a  role,  it  has 
not  been  established  under  what  conditions  certain  effects  need  to 
be  included  in  analyses.  The  most  important  of  these  appear  to  be 
contact  deformation  between  impactor  and  structure,  transverse 
shear  deformations,  viscoelastic  material  behavior,  and  large 
deformation  (membrane)  effects.  The  research  reported  here  was 
undertaken  to  determine  the  relative  importance  of  these  phenomena 
in  predicting  the  impact  response  of  laminated  composites. 

There  are  several  types  of  blunt  object  impact  damage  which  can 
occur  in  laminated  composites  [1-3]*: 

-  penetration  caused  by  high-velocity  impact; 

-  crushing  of  impacted  surface  material,  delamination,  and 
shattering  of  the  back  surface  at  moderate  velocities  which  is 
caused  by  the  interaction  of  elastic  wave  effects  and 
quasistatic  structural  stresses;  and 

-  a  combination  of  flexural  cracking  and  interply  damage  which 
can  occur  at  low  impactor  velocities  and  may  not  involve  wave 
propagation  effects. 

It  is  the  incipient  damage  caused  by  low  velocity,  blunt  object, 
transverse  impact  with  which  this  report  is  concerned. 


*Numbers  in  square  brackets  identify  references. 
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Several  investigators  have  researched  the  problem.  Chou  and 
Flis  [4]  have  shewn  that  analysis  of  beam  and  plate  structures 
considering  only  flexure  and  static  deflection  shapes  can  give 
non-conservative  predictions  of  structural  strains.  They  also  show 
[4,  5]  that  solving  problems  dynamically  but  still  treating  only 
flexure  as  was  done  by  McQuillen,  Llorens,  and  Gause  [6]  and  Hayes 
and  Rybicki  [7]  significantly  overpredicts  the  strains.  Llorens, 
McQuillen,  and  Gause  [8,  9,  10]  attempted  to  correct  this 
overprediction  using  exact  and  approximate  viscoelastic  damping, 
but  were  only  partially  sucessful:  non-conservative  results  were 
obtained  in  several  cases.  It  was  demonstrated  [4]  that  treatment 
of  Herzian  contact  deformation  between  impactor  and  structure  as 
was  done  by  Sun  and  Chattopadhyay  [11],  Bostaph  and  Elber  [12],  and 
Elber  [13]  also  tends  to  correct  the  elastic  flexural 
overpredictions,  but  direct  'comparison  with  experimental  data  was 
not  available. 

There  is  an  uncertainty  whether  viscoelastic  structural  damping 
or  Hertzian  contact  deformation  is  more  important  to  achieving 
accurate  predictions  of  structural  strains  in  low-velocity  impact 
of  composite  beams  and  plates.  References  [12]  and  [13]  show  that 
membrane  force  effects  can  also  play  an  important  role  in  obtaining 
accurate  predictions  of  stresses  and  deformations  in  thin  composite 
plates  undergoing  impact.  It  is  also  well-known  that  transverse 
shear  deformations  can  be  important  for  thick  plates.  Several 
researchers  [10,  14-16  for  example]  have  used  finite  element 
methods  to  analyze  composite  laminate  impact  problems,  some 
including  the  effects  of  elastic  contact  stress  [14,  15]  ,  contact 
deformation  [16]  and  transverse  shear  deformation  [14-16].  Several 
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of  the  methods  were  semi-empirical  or  approximate  and  have  not  been 
substantiated  by  experimental  evidence.  None  has  investigated  the 
relative  importance  of  structural  viscoelasticity,  contact 
stresses,  transverse  shear  deformations,  and  membrane  effects. 

A  program  of  analysis  and  experiment  was  initiated  to 
investigate  under  what  conditions  each  of  the  four  phenomena 
described  above  may  be  important  to  analyzing  composite  structures. 
The  program  consisted  of: 

1.  An  investigation  into  viscoelastic  properties  of  composites 
to  determine  which  were  the  most  likely  modes  of  viscoelastic 
dissipative  mechanisms  during  impact  and  to  determine  the 
short-time  viscoelastic  properties  of  AS-3501  graphite /epoxy. 

2.  Finite  element  analysis  of  low-velocity  impact  of  composite 
plates . 

3.  Impact  tests  of  AS-3501  graphite/ epoxy  plates. 

The  original  plan  was  to  analyze  AS  3501  graphite /epoxy 
composite  beams  using  the  FEAP74  finite  element  code  [17-23].  The 
code  was  developed  for  the  U.  S.  Department  of  Transportation  for 
impact  studies  and  contains  3-D,  2-D,  shell  (with  transverse 
shear),  and  contact  elements;  viscoelasticity;  orthotropy;  and 
large  displacements.  When  received,  the  program  was  found  to  be 
defective  in  the  viscoelasticity  and  large  deflection  shell 
elements.  Also,  the  viscoelastic  element  was  not  orthotropic.  As 
a  result,  the  analysis  portion  of  the  program  was  limited  to 
elastic  finite  element  analysis  which  treated  elastic  shear  and 
contact  deformations.  The  analyses  were  used  to  evaluate  the 
relative  importance  of  these  effects  on  prediction  of  strains  and 
to  provide  an  analytical  baseline  for  comparison  with  tests. 
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In  order  to  evaluate  the  accuracy  of  elastic  analysis 
predictions  for  behavior  of  laminated  composites  undergoing  impact, 
tests  were  conducted  on  specially  fabricated  AS-3501  graphite/epoxy 
plates  using  the  NAVAIRDEVCEN  instrumented  impact  tower.  Results 
of  the  tests  were  compared  with  the  finite  element  analyses,  and 
conclusions  were  drawn  concerning  the  need  for  including  transverse 
shear,  contact,  and  other  deformation  phenomena  in  an  impact 
analysis  capability. 

This  report  presents  analysis  and  test  work  performed  under  the 
program,  recommends  a  viscoelastic  model  of  AS-3501  composite 
plates  for  use  in  impact  analysis,  and  draws  conclusions  concerning 
the  importance  of  elastic  and  viscoelastic  behavior  to  impact 
analysis  of  composite  laminates. 
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II.  VISCOELASTIC  BEHAVIOR  OF  COMPOSITE  LAMINATES 

A.  DISSIPATIVE  MECHANISMS  DURING  IMPACT 

The  work  of  Llorens,  McQuillen,  and  Gause  [3  -  10]  has 
demonstrated  that  an  energy  absorption  mechanism  other  than  elastic 
flexure  is  active  during  low-velocity  impact  (less  than  9  m/s  or  30 
ft/s).  Their  analysis  and  comparison  with  test  data  for  impacted 
AS-3501  plates  show  that  flexural  viscoelastic  analysis  compares 
much  better  with  experimental  results  than  purely  elastic  flexural 
analysis.  However,  there  are  still  large  discrepancies  for  certain 
plate  aspect  ratios  indicating  that  the  energy  absorption  mechanism 
may  not  be  flexural,  and  that  perhaps  a  different  viscoelastic 
model  might  improve  accuracy  of  analytical  predictions. 

1.  Flexural  Viscoelasticity  of  AS-3501  Graphite/Epoxy 

It  is  well-known  [24-27,  for  example]  that,  for  a 
unidirectional  composite  consisting  of  high-stiffness  elastic 
fibers  in  a  viscoelastic  polymeric  matrix,  the  effective 
viscoelastic  composite  behavior  is  most  pronounced  under  transverse 
normal  stress  022'  transverse  shear  stress  T ^3  andT23,  and  axial 
shear  stress  T 12  (Figure  1).  Viscoelastic  behavior  will  be 
exhibited  under  axial  normal  stress,  but,  due  to  the  dominance  of 
the  high-stiffness  elastic  fibers,  it  can  be  negligible  compared  to 
the  other  modes.  In  plate  bending,  however,  the  mode  of  deformation 
can  be  predominately  flexure  and  it  may  be  that  axial  normal 
viscoelasticity  alone  can  be  great' enough  to  explain  test  results 


on  impacted  plates. 
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In  order  to  see  if  viscoelastic  behavior  of  AS-35G1 
graphite/ epoxy  under  axial  normal  stress  is  great  enough  to  give 
the  flexural  viscoelastic  behavior  required  to  correlate  analysis 
with  test  data  [8  -  10],  viscoelastic  properties  of  AS-3501 
obtained  from  the  tests  of  Renton  and  Ho  [28]  were  analyzed  and 
compared  with  the  flexural  viscoelastic  properties  used  in 
references  [8  -  10]. 

a.  Matrix  Shear  Power  Law  Creep  Compliance  [28].  -  Renton  and 
Ho  [28]  have  done  tensile  creep  tests  of  [±45]s  laminates  of 
AS-3501  and  have  used  a  power  law  creep  compliance  given  by 

e/s  =  D0  +  D^tn  (1) 

where 

e  =  tensile  strain 
s  =  tensile  stress 
t  =  time 

Dq,  D^,  and  n  =  viscoelastic  constants. 

Under  all  environments,  n  =  0.18  gave  good  correlation  with 
test  data.  The  constants  and  Dq  for  room  temperature  tests  at 
50%  relative  humidity  were  found  to  be 


D0  =  3 . 412 ( 10“7 )  in2/lb  (4.949(10“5)  MPa-1} 

Dx  =  0 . 108 ( 10"7 )  in2/lb  {0.155(10"5)  MPa"1} 

Using  laminate  analysis  and  assuming  that  the  Halpin-Tsai 
equations  relating  constituent  properties  to  composite  properties 
are  valid,  Renton  and  Ho  construct  a  nomogram  by  which  the  matrix 
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shear  creep  compliance,  FM,  can  be  obtained  from  [±45]  tensile 
creep  compliance.  Using  the  values  of  Dg  and  D]_  given  above  and 
Figure  24  of  Reference  [28],  the  following  relationship  is  obtained 
for  room-temperature  matrix  shear  creep  compliance: 

Fm  =  9.30{10"6)  +  0 . 340 ( 10“6 ) t°  * 18  in2/lb  (2) 

=  13.5(10-4)  +  0 . 493 ( 10’4)t0*18  MPa-1 

b.  Matrix  Extensional  Power  Law  Creep  Compliance.  -  Even  if 
matrix  isotropy  is  assumed,  one  must  know  the  bulk  creep  compliance 
or  the  creep  Poisson  compliance  to  be  able  to  determine  the 
extensional  creep  compliance,  Dj^  of  the  matrix  [27].  However,  for 
an  order-of -magnitude  estimate,  it  is  assumed  that  Poisson's  ratio 
vM  remains  constant  with  time  at  a  value  of  0.3.  For  isotropic 
behavior,  this  results  in  an  extensional  creep  compliance  for  the 
matrix  of 

Dm  =  0.5Ff,j/(l  +  Vjkj) 

=  3 . 58 ( 10“6 ) ( 1  +  0. 03656t° • 18 )  in2/lb  (3) 

=  5 . 19 ( 10~4 ) ( 1  +  0.03656t0-18)  MPa’1 
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For  AS-3501,  fiber  volume  fraction  is  VF  =  0.62.  Fiber  modulus  can 
be  calculated  from  equation  (3)  with  t  =  0  and  equation  (4)  with 
composite  initial  modulus  1/D*  =  21.0  Mlb/in2  (145  GPa).  The 
result  is  1/Dp  =  33.75  Mlb/in2  (232.7  GPa).  For  arbitrary  time, 

( 4 )  gives 

% 

D*  =  4.75(10-8)(1  +  0 . 03656t° • 18  )  / ( 1  +  0 . 03637t° • 18 )  in2/lb  (5) 
=  6.90(10-6)(1  +  0 . 03656t° • 18 ) / ( 1  +  0 . 03637t° • 18 )  MPa-1 

as  the  axial  creep  compliance  of  unidirectional  AS-3501. 

d.  Laminate  Extensional  Creep  Compliance  Power  Law 
Approximation.  -  In  order  to  compare  the  viscoelastic  properties 
of  AS-3501  to  those  used  for  impacted  plate  experimental 
correlation  in  [8-10],  it  is  noted  that  the  extensional  moduli  of 
the  multidirectional  laminates  used  in  [8-10]  are  approximately 
half  that  of  the  unidirectional  composite.  Assuming  that  equation 
(4)  will  approximate  the  shear  creep  compliance  of  the 
multidirectional  plate  if  2DF  is  used  in  place  of  Dp,  the  following 
equation  is  obtained  for  an  order-of-  magnitude  estimate  of  plate 
extensional  creep  compliance  Dp*: 

Dp*  =  9 . 46 ( 10-8 ) ( 1  +  0.03656t°-ls)/(l  +  0 . 03619t° • 18 )  in2/lb  (6) 
=  13.72U0-6)(1  +  0.03656t°-18)/(l  +  0 . 03619t° • 18 )  MPa-1 

Another  order-of-magnitude  estimate  would  be  to  double  the 
unidirectional  axial  creep  compliance  (5).  This  results  in  an 
approximation  which  is  numerically  very  close  to  equation  (6); 
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therefore,  (6)  will  be  used  as  the  laminated  plate  extensional 
creep  compliance. 

e .  Exponential  Approximation  of  Laminate  Extensional  Creep 
Compliance .  -  Since  the  viscoelastic  creep  function  assumed  in  the 
plate  analysis  of  references  [3-10]  is  an  exponential  law,  it  is 
desirable  to  take  the  power  law  creep  compliance  (6)  obtained  from 
Renton  and  Ho's  data  and  fit  it  to  an  exponential  three-parameter 
viscoelastic  solid  having  a  creep  relaxation  function  of 


dP*RH  =  1/<3l  +  (P/92  "  l/qi)exp[-g1t/q2]  (7) 

where  q^,  q2,  and  p  are  viscoelastic  constants.  Fitting  (7)  to  (6) 
at  t  =  0  s,  1.0s,  and  infinity  gives  the  following  values  for  the 
constants 

q1  =  10.46  Mlb/in2  (72.13  GPa) 

q2  =  294.  Mlb.s/in2  (2029.  GPa.s)  (8a,b,c) 

p  =  27.85  s 


It  is  noted  that  these  creep  compliance  calculations  are 
approximate  and  are  to  be  used  for  order-of-magnitude  comparisons 
only! 


2 .  Flexural  Viscoelasticity  Required  to  Correlate  Impact  Data 


The  flexural  viscoelastic  constitutive  relation  usee,  by 


Llorens,  McQuillen,  and  Gause[8-10]  to  correlate  plate  impact 
data  is  equivalent  to  a  Kelvin  solid  having  the  following 


extensional  creep  compliance: 
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°P*LMG  =  (l/q3)(l  -  exp[ -q3t/q4])  (9) 

which  has  the  form  of  (7),  but  with  p  =  0.  Their  analysis  of 
impacted  plates  correlates  with  tests  of  AS-3501  laminates  when 


q4/q3  =  25  { 10-6  )  s  to  100(10_6)s. 


(10a) 


For  an  accurate  comparison  between  (9)  and  (7),  is  necessary  that 


q3  =  qi  (10b) 

In  order  to  compare  the  amount  of  extensional  viscoelasticity 
required  by  Llorens,  McQuillen,  and  Gause  to  correlate  their  test 
data  with  the  extensional  viscoelasticity  available  from  AS-3501  as 
determined  above  from  the  tests  of  Renton  and  Ho,  it  is  desirable 
to  have  a  quantitative  measure  of  damping  which  is  not  heavily 
dependent  upon  the  form  of  the  assumed  viscoelastic  constitutive 
relation.  The  amount  of  energy  per  unit  volume  dissipated  in  one 
cycle  of  sinusoidal  loading,  EDC,  is  one  such  convenient  quantity 
which  is  calculated  from  the  equation  [29,  30] 

EDC  =  T  w(qB  -  qAP)So2/(qA2  +  <?B2w2)  (11) 

where 

w  =  loading  frequency 

q^,  qB,  p  =  viscoelastic  constants  ( A  =  1  or  3 ,  B  =  2  or  4 

sq  =  stress  half -amplitude 


The  ratio  of  energy  density  dissipated  per  cycle  using  the  material 
behavior  required  to  correlate  plate  impact  data  (EDCLMG)  to  that 
which  has  been  calculated  using  known  AS-3501  material  behavior 
( EDCt^tj )  becomes 

q4(qi2  +  q22w2) 

EDClmg/EDCrh  =  - - - — - "  (12) 

(q3^  +  q^w*)  (q2  -  q]_p) 

and  should  tell  if  flexural  viscoelasicity  or  some  other  energy 
absorption  mechanism  is  dominating  impact  behavior. 

The  tests  of  Chou,  Flis,  and  Miller  [31]  were  used  to  correlate 
the  analysis  of  [8-10].  These  tests  had  typical  impact  event  times 
of  5  ms  which  are  equivalent  to  an  oscillatory  period  of  10  ms  and 
a  cyclic  frequency  w  of  approximately  600  rad/s.  Using  this 
frequency  and  a  representative  value  of  q4/q3  =  60  s,  equations 
(8),  (10),  and  (12)  give 

EDCljijg/EDCrh  =  6(10^)  (13) 

Clearly,  the  degree  of  flexural  viscoelastic  behavior  required 
to  fit  impact  data  is  several  orders  of  magnitude  greater  than  that 
which  exists  in  the  material.  It  is  concluded  that  the  energy 
absorbing  mechanism  cannot  be  flexural  viscoelasticity.  The 
determination  and  modeling  of  the  true  cause(s)  of  the  disparity 
between  elastic  impact  analysis  and  test  data  may  provide  better 
correlation  between  analysis  methods  and  tests  and  yield  guidelines 
for  more  accurate  design  approaches  to  impact  resistance. 
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3 .  Impact  Snercry  Absorbing  Mechanisms 

Energy  absorbed  by  composite  plates  during  impact  is  not 
entirely  due  to  elastic  flexure,  and  Llorens,  McQuillen,  and  Gause 
have  shown  that  other  mechanisms  are  operative.  The  preceeding 
analysis  has  shown  that  flexural  viscoelasticity  is  negligible  and 
that  other  energy  absorptive  phenomena  must  be  treated  in  order  to 
obtain  good  analytical  predictions  for  design  purposes.  As 
presented  in  the  introduction  to  this  report,  the  following 
mechanisms  are  considered  to  be  the  strongest  candidates: 

1.  Contact  deformations  between  impactor  and  plate. 

2.  Transverse  shear  deformations. 

3.  Membrane  forces. 

Of  these,  the  membrane  force  effects  are  known  to  be  important 
whenever  plate  deflections  exceed  plate  thickness  (methods  of 
treatment  of  large  deflections  of  thin  composite  laminates  during 
impact  analysis  have  been  described  by  Elber  and  Bostaph  [12,13], 
for  example).  Shorter,  thicker  plates  appear  to  require  treatment 
of  contact  deformations  and/or  transverse  shear  deformations. 

Both  contact  deformations  and  transverse  shear  deformations  are 
matrix-dominated  phenomena,  i.e.,  the  properties  of  the  matrix 
material  will  control  laminate  behavior  in  these  modes. 
Fiber-dominated  laminates  are  therefore  stiff  in  flexure  (or  any 
other  in-plane  deformation  mode  such  as  membrane  stretching) ,  but 
may  be  relatively  flexible  under  transverse  normal  stress  and 
transverse  shear  stress.  Transverse  normal  stress  effects  will 
control  contact  deformations  induced  by  hard  impactors.  For 
composites,  transverse  shear  effects  can  be  important  for 
relatively  thin  plates.  For  example,  it  is  well-known  that  the 
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deflection,  u,  of  a  cantilevered  wide  plate  under  a  concentrated 
transverse  load,  P,  at  its  tip  can  be  approximated  by 
strength-of-materials  analysis  as 

u  =  ( PL3 / 3CI ) ( 1  +  0.25Ch2/GL2) 

where 

L  =  plate  length 
h  =  plate  thickness 
I  =  wh3/12,  w  =  plate  width 
G  =  material  transverse  shear  modulus 
C  =  plate  extensional  modulus 

The  last  term  in  the  last  parentheses  represents  the  contribution 
of  transverse  shear  deformation  to  deflection.  For  isotropic 
materials  where  E  and  G  are  the  same  order  of  magnitude,  the 
transverse  shear  deflection  will  be  no  more  than  3%  of  the  total  if 
length  L  is  longer  than  5h.  For  a  typical  graphite  epoxy,  L  must 
be  greater  than  2 Oh  for  the  same  relative  magnitude  of  shear 
deflection. 

In  addition,  since  polymeric  matrix  behavior  is  viscoelastic, 
it  is  possible  that  treatment  of  elastic  transverse  normal  and 
shear  is  insufficient:  viscoelastic  contact  deformations  and 
transverse  shear  deformations  may  need  to  be  modeled  in  order  to 
obtain  reasonable  analytical  accuracy  for  design  or  analysis 


purposes . 
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B.  VISCOELASTIC  PROPERTIES  OF  AS- 3 5 01 


1.  Viscoelastic  Material  Modelinc 


It  has  been  demonstrated  here  and  elsewhere  that  in-plane 
stress-strain  relations  can  be  elastic  without  loss  of  accuracy  for 
low- velocity  plate  impact  analysis.  In  order  that  a  plate  exhibit 
contact  deformation  viscoelasticity  and  transverse  shear 
viscoelasticity,  it  will  be  necessary  to  model  transverse  normal 
and  shear  stress-strain  relations  as  viscoelastic.  An  elementary 
model  which  would  take  the  major  effects  into  account  although  it 
may  not  be  thermodynamically  correct  is  to  assume  that  composite 
transverse  normal  and  shear  viscoelastic  relations  are  uncoupled 
and  have  the  following  exponential  relaxation  functions: 


<%>  ~  G0  +  G^expf-t/b)  (shear) 
CT  =  C0  +  C^expf-t/a)  (normal) 


(14) 


where  Cq,  Gq,  C^,  G]_,  a,  and  b  are  constants. 

This  would  allow  use  of  certain  finite  element  codes  which  use  this 
functional  form  ([17],  for  example,  with  appropriate  code 
modification  for  orthotropy) . 

In  order  to  obtain  the  three  material  constants  required  in 
each  of  equations  (14)  and  (15),  tests  must  be  performed.  The  most 
efficient  way  to  obtain  shear  properties  for  equation  (14)  is  to 
perform  a  creep  test  of  a  [±45JS  laminate  in  tension  as  was  done  by 
Renton  and  Ho  [28].  The  creep  compliance  thus  obtained  can  be 
inverted  to  give  the  relaxation  modulus  (14)  by  Laplace  transform 
techniques  using  the  correspondence  principle  of  linear 


-15- 


MADC-37'' 06-60 


viscoelasticity  [27,  for  example].  Using  relationships  developed 
between  matrix,  fiber,  and  composite  viscoelastic  stress-strain 
relations  [24,  27],  one  may  determine  matrix  viscoelastic 
stress-strain  relations  from  composite  shear  behavior  and  in  turn 
use  the  matrix  behavior  to  calculate  the  transverse  normal 
viscoelastic  relaxation  modulus  ! 15 )( see  reference  [28]  for  one 
description  of  this  process). 

Experimentally,  then,  creep  tests  on  [±45]s  tension  specimens 
are  sufficient  to  determine  any  viscoelastic  properties  that  may  be 
important  to  impact  analysis  of  composite  laminates.  All 
multidirectional  laminates  of  the  same  unidirectional  material  have 
approximately  the  same  transverse  normal  and  shear  behavior  since 
stacking  sequence  affects  only  in-plane  constitutive  relations. 

The  time  from  impact  to  rebound  or  fracture  will  depend  upon 
material  stiffness  properties,  structural  geometry,  and  impactor 
mass.  During  low-velocity  impact  of  aircraft  skins  by  stones, 
tools,  or  other  similar-size  masses,  this  elapsed  time  is  seldom 
greater  than  0.05  s  and  can  be  considerably  shorter.  The 
viscoelastic  test  data  for  [±45}  AS-3501  graphite /epoxy  laminates 
generated  by  Renton  and  Ho  [28]  used  creep  data  intervals  of  one 
minute  (60  s).  Since  the  impact  events  are  taking  three  orders  of 
magnitude  less  time  than  the  first  data  point  in  [28],  the  use  of 
Renton  and  Ho’s  data  for  accurate  viscoelastic  deformation 
calculations  is  questionable. 

The  following  section  describes  tests  performed  to  determine 
the  shear  relaxation  modulus  of  AS-3501  graphite/epoxy  for  use 
during  short  time  intervals  typical  of  impact  occurrences. 
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2 .  Test  Program 

In  order  to  develop  viscoelastic  properties  which  might  be 
valid  over  time  intervals  on  the  order  of  fractions  of  a  second,  a 
series  of  high-speed  creep  tests  was  performed  on  [-45]  23  AS-3501 
using  an  Instron  servohydraulic  tensile  test  machine  at  the 
laboratories  of  the  Naval  Air  Development  Center. 

Each  specimen  was  0.5  in  (25.4mm)  wide  and  9.0  in  (225  mm) 
long  and  was  appropriately  tabbed  with  glass  end  tabs.  Each  was 
instrumented  with  a  3-element  rectangular  strain  gage  rosette  with 
the  middle  gage  in  the  direction  of  loading  and  the  other  two  at 
angles  of  +45°  and  -45°  to  the  loading  direction. 

Load  and  axial  strain  were  read  into  a  Nicolet  digital  storage 
scope,  and  the  45°  strain  gages  were  monitored  by  a  Tektronix 
fluorescent  storage  oscilloscope.  Photographic  records  of  45 
strains  were  made  with  a  Polaroid  camera  fitted  to  the  Tektronix 
scope . 

Specimens  were  "instantaneously"  ramp- loaded  to  a  predetermined 
load  (ranging  from  125  lb  [556  N]  to  640  lb  [2.85  kN])  and  held  at 
this  constant  load  for  approximately  ten  minutes.  During  the  first 
twenty  seconds,  2,000  data  points  were  recorded  by  the  Nicolet 
storage  scope  for  both  load  and  strain.  This  provided  10  data 
points  in  the  first  0.1  s  which  is  considerably  better  than  one 
data  point  every  60  s.  Twenty-second  data  records  were  also  made 
after  approximately  five  minutes  and  ten  minutes,  respectively,  to 
obtain  long-time  creep  information. 
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3 .  Test  Results 

Typical  creep  test  data  are  presented  in  Figures  2,  3,  and  4 
which  give  axial  load  and  strain  for  the  first  twenty  seconds,  the 
twenty  second  interval  beginning  after  an  elapsed  time  of  five 
minutes,  and  the  twenty  second  interval  after  an  elapsed  time  of 
thirteen  minutes,  respectively. 

It  was  found  that  during  the  initial  loading,  approximately  two 
data  intervals  of  0.01  s  were  required  for  the  test  machine  to 
reach  the  desired  load,  and  another  28  data  intervals  for  the 
hydraulics  to  stabilize  the  load  (one  or  two  "dips"  of  about  4% 
from  the  preset  load  typically  occurred) .  The  strain  data  was 
adjusted  for  the  load  variations  by  multiplying  the  strain  at  a 
given  time  by  the  ratio  of  actual  load  at  that  time  to  the  stable 
load  reached  after  0.5  s.  Also,  as  is  evident  in  Figure  2,  noise 
was  superimposed  on  the  strain  signal.  Data  were  assumed  to  lie 
midway  between  noise  peaks.  Figure  5  shows  the  results  of 
adjusting  and  smoothing  the  data  from  Figure  2  (note  the  expanded 
time  scale  in  Figure  5 ) . 

Adjusted  extensional  strain  data  in  the  0°,  +45°,  and  -45° 
directions  from  the  strain  gage  rosette  were  transformed  using  the 
common  plane  transformation  equations  to  yield  the  axial  shear 
strain  (parallel  to  fiber  directions).  Load  was  divided  by  the 
[±45]  specimen  cross-sectional  area  to  obtain  tensile  stress  which 
was  then  transformed  to  axial  shear  stress  (parallel  to  fiber 
direction  -  Figure  1).  The  resulting  shear  stress  and  time  history 
of  shear  strain  were  analyzed  to  determine  axial  shear  relaxation 
modulus  as  follows: 


W, 

& 


Figure  4.  Load  and  axial  strain  vs.  time  dur 
[ ±45 ]  AS-3501  graphite /epoxy  tensL 
lb  (543  N)  load.  Time  interval  78 
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Figure  5.  Adjusted  axial  strain  vs.  time  for  creep  test  of  [±45]  AS-3501 

graphite/epoxy  tension  specimen  under  125  lb  (556  N)  load.  Time 
interval  0.0  s  to  20  s. 
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The  form  of  shear  relaxation  modulus  chosen  was  that  of 
equation  (14)  which  has  three  viscoelastic  constants  -  the  static 
(slow  load)  shear  modulus  Gq,  the  viscous  modulus  G]_,  and  the 
exponential  decay  parameter,  b: 


g12  =  G0  +  Giexp(-t/b)  (16) 

It  can  be  shown  that  the  shear  creep  compliance  corresponding  to 
this  relaxation  modulus  is 


F12  =  1  /Gq  -  G^Ggt  Gq  +  G]_ )  ] -1exp  ( -t/c  ) 


where 


C  —  t(Gg  +  Gj_)/GQ 


(17a) 

(17b) 


The  static  shear  modulus  Gq  was  determined  from  slow-speed  tension 
tests  of  specimens  similar  to  those  used  for  creep  tests.  The 
ratio  of  viscous  modulus  G^  to  shear  modulus  Gq  was  calculated  from 
the  equation 


Gi/G Q  =  Gjj/Gq  -  1 


(18) 


where  G^  is  the  dynamic  shear  modulus  found  from  high-speed  test 
data.  In  this  case,  G^  was  determined  from  plots  of  stress  versus 
strain  obtained  from  the  initial  0.05  s  of  creep  data.  With  Gq  and 
G^  known,  the  exponential  creep  decay  parameter,  c,  was  found  by 
choosing  one  point  on  the  adjusted  extensional  strain  curve  and 
solving  for  c  from  equation  (17a).  The  exponential  relaxation 
decay  parameter  b  was  then  determined  from  equation  (17b). 
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The  parameter  b  was  calculated  from  several  points  on  the 
initial  portion  of  the  creep  strain  curve  and  was  found  to  be 
affected  by  the  point  chosen  for  its  determination.  The  most 
consistent  results  were  obtained  when  data  points  were  taken  within 
the  the  first  0.3  s  of  creep  data  where  the  majority  of  the  creep 
strain  took  place. 

Computed  relaxation  modulus  constants  obtained  from  the 
viscoelastic  creep  tests  are  presented  in  Table  1.  Appendix  A 
contains  computer-generated  plots  of  the  creep  data  from  the  tests. 


Table  1.  Relaxation  Modulus  Constants  for  AS-3501  in  Axial  Shear. 


Relaxation 

Modulus 

Constant 

Value 

[  Gi  2  =  Go  +  G]exp(-t/b)] 

G0 

0.941(10®)  lb/in2  [6.49  MPa] 

Gi/Gq 

0.0313 

b 

0.15  +  0.05  s*1 

I 
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III.  ELASTIC  ANALYSIS  OF  IMPACTED  COMPOSITE  BEAMS 


A.  FINITE  ELEMENT  MODEL 


1 .  Structural  Configuration  and  Material  Properties 

Elastic  analysis  of  AS-3501  graphite/epoxy  composite  plate 
structures  was  performed  to  determine  the  relative  effects  of 
contact  deformations  and  transverse  shear  deformations,  and  to 
provide  an  analytical  baseline  for  comparison  of  data  from  impact 
tests.  For  both  analytical  and  experimental  simplicity,  a  long 
plate  clamped  at  both  ends  was  chosen.  Plates  ranged  in  length 
from  2.56  in  (65  mm)  to  13.98  in  (355  mm),  but  all  plates  were  1.57 
in  (40  mm)  wide.  Plates  had  thicknesses  of  1/8  in  (3  mm)  or  1/4  in 
(6  mm)  corresponding  to  24-  and  48-plies  of  AS-3501,  respectively. 
Figure  6  illustrates  plate  configurations  analyzed. 

The  actual  stacking  sequences  used  in  impact  tests  wv:ich  the 
finite  element  analysis  was  designed  to  model  were 
[ (±45/02)2/±45/0/90]s  and  [ (±45/02)2/±45/0/90]2s.  For  the  model, 
each  plate  was  assumed  to  have  uniform  orthotropic  elastic 
constants  corresponding  to  [±45/02 ^NS  stacking  sequences  of  AS-3501 
graphite/ epoxy.  Laminate  analysis  confirmed  that  axial  strains  in 
a  plate  contained  less  than  1%  error  with  the  uniform  material 
assumption  compared  to  individual  ply  modeling. 

The  major  material  axis  was  chosen  to  be  oriented  at  either  0° 
or  90°  to  plate  longitudinal  axis  for  each  of  the  four  plate 
geometries,  which  gave  a  total  of  eight  plates  to  be  analyzed. 
Figure  6  presents  geometry,  material,  and  associated  numbering 
system  details. 


WIDTH  W  THKNS  h 
in/mm  in/mm 


STACKING  SEQUENCE 
(X  =  0°) 


MATERIAL 

CODE 


14/355 


5.31/135 


2.56/65 

»1 


t  (±45/02)2/±45/0/90]s 
[ (±45/902T2/±45/90/0Ts 
[ {±45/02T2/±45/0/9 0]s 

[ (±45/902T2/±45/90/0Ts 

[ (±45/02T27±45/0/90]>s 
[ (±45/902)2/±45/90/0J2s 
[ (±45/02)2/±45/0/90]2s 
[ (±45/902)2/±45/90/oT2S 


>/±45/90/0]s 
'±45/0/9 0]s 
)/±45/90/0Ts 


MATERIAL  CODE 

AVERAGE  PLANE  STRAIN  ELASTIC  CONSTANTS* 

Mpsi  (GPa) 

CXX 

CZZ 

CXZ 

CSS 

H-  High  Stiffness 
in  X  Direction 

11.2(75.4) 

1.71(11.8) 

0.86(5.91) 

1.30(8.96) 

L-  Low  stiffness 
in  X  Direction 

5.45(37.6) 

1.71(11.8) 

.712(4.91) 

1.30(8.96) 

*  ^Normal  Stress  s  ®XX  —  ^XX  ®XX  ^XZ  ®ZZ 

I  and  Strain  e: 

szz  =  cxz  exx  +  czz  ezz 

Shear  Stress  s^z  and  Strain  eyz:  sxz  =  CSS  eXZ 
(Note  that  evZ  is  tensorial  shear  strain) 


Figure  6.  Structural  and  material  configurations  analyzed  by 
FEAP74  finite  element  code. 
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Since  plates  were  an  order  of  magnitude  wider  than  they  were 
thick,  a  two-dimensional  plane  strain  analysis  was  appropriately 
selected  for  the  impact  analyses.  In  order  to  determine  the 
effects  of  elastic  contact  deformations  and  transverse  shear 
deformations  on  analysis  results,  transverse  material  properties 
were  altered  and  finite  element  analyses  were  run  in  the  following 
sequence : 

a.  full  two-dimensional  orthotropic  elasticity  which  contained 
the  effects  of  flexure,  elastic  contact  deformations,  and  elastic 
transverse  shear  deformations. 

b.  transverse  normal  stiffness  was  increased  by  two  orders  of 
magnitude  which  eliminated  contact  deformations  but  maintained 
flexure  and  transverse  shear  effects. 

c.  both  transverse  normal  and  transverse  shear  stiffnesses 
were  incraesed  by  two  orders  of  magnitude  which  eliminated  contact 
and  transverse  shear  deformations  but  maintained  flexural  effects. 

2 .  Finite  Element  Code 

The  finite  element  code  chosen  for  the  analysis  was  FEAP74,  a 
program  developed  by  University  of  California,  Berkeley,  for  the  U. 
S.  Department  of  Transportation  especially  for  contact-impact 
problem  analysis.  The  code  contains  operational  two-  and 
three-dimensional  orthotropic  elastic  elements,  laminated 
orthotropic  plate  element,  and  contact  elements.  Dynamic  impact 
problems  can  be  analyzed  using  explicit  or  implicit  finite 
difference  formulations  in  the  time  domain.  The  implicit  option 
was  chosen  since  results  thus  obtained  are  always  stable  and 
convergent. 
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Unfortunately,  large  deflection  plate  and  viscoelastic  elements 
were  inoperable  in  the  version  received  from  Cal-Berkeley .  while 
these  elements  were  not  able  to  be  made  operational,  several  other 
errors  were  found  in  the  code  and  corrected  so  that  the  other 
elements  will  run  with  the  implicit  dynamic  analysis  option.  The 
version  currently  on  tape  storage  at  the  NADC  computer  facility 
contains  all  corrections  which  were  made. 

3 .  Finite  Element  Mesh  Construction 

In  order  to  prevent  possible  numerical  instabilities  during  the 
finite  element  analyses,  care  was  taken  to  construct  elements 
having  nearly  equal  stiffnesses  in  the  principal  directions.  Plane 
strain  meshes  sized  from  64  elements  and  87  nodes  to  383  elements 
and  462  nodes  were  analyzed  for  optimal  accuracy  and  running  time. 
It  was  found  that  less  than  3%  difference  in  pertinent  stresses  and 
strains  was  obtained  between  meshes  having  four  elements  through 
the  beam  thickness  and  six  elements  through  the  beam  thickness. 
Accordingly,  four-element-thick  beam  meshes  were  used  for  the 
finite  element  production  runs. 

The  impactor  from  the  NADC  impact  tower  required  less  precision 
in  modeling  than  did  the  plates.  Force  in  the  load  cell  of  the 
impactor  was  the  quantity  of  interest  since  time  histories  of  the 
load  could  be  obtained  during  tests.  Care  was  taken  to  model  the 
impactor  tip  so  that  accurate  contact  stresses  would  be  obtained, 
but  the  remaining  parts  (load  cell,  tip  connector,  frame,  and 
guides)  were  modeled  as  axial  elements.  A  dynamic  study  of  the 
impactor  frame  was  made  to  determine  how  accurate  a  dynamic  model 
of  the  frame  was  necessary.  Strength-of-materials  analysis  showed 
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that  natural  frequencies  of  the  frame  and  its  parts  could  be  on  the 
same  order  of  magnitude  as  those  of  the  plates,  and  that  careful 
attention  needed  to  be  paid  to  stiffnesses  as  well  cs  masses  of  the 
frame  parts.  The  resulting  impactor  mesh  was  therefore  designed  to 
be  dynamically  representative  during  impact. 

Figure  7  illustrates  the  mesh  developed  for  the  17  beam  and  the 
8.4  lb  (3.8  kg)  impactor  which  has  a  1/8-in-  (3.2-mm-)  radius 
cylindrical  tip.  The  remaining  beam  meshes  are  similar. 

B.  RESULTS 

Finite  element  analysis  results  were  obtained  for  ail  eight 
plate  configurations  with  (a)  full  elastic  behavior,  (b)  flexure 
and  transverse  shear  deformations  only,  and  (c)  flexure  only.  The 
longest  plates,  II  and  12,  with  aspect  ratios  ( length-to-thickness ) 
of  about  100,  suffered  numerical  instability  problems  when  both 
transverse  shear  and  normal  stiffnesses  were  artificially  increased 
together.  An  attempt  at  resolving  the  problem  by  creating  a  finer 
mesh  was  thwarted  by  size  limitations  imposed  by  the  FEAP74  code. 
All  other  runs  were  completed  and  results  are  available. 

Runs  were  made  at  impactor  initial  velocities  ranging  from  42 
in/s  (500  mm/s)  to  315  in/s  (8  m/s).  It  was  determined  that 
impactor  force  and  beam  flexural  strain  results  were  proportional 
to  velocity,  and  therefore  only  the  315  in/s  (8  m/s)  results  are 
presented  here. 
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Figure  7.  Finite  element  mesh  for  2,56  in  (65  mm)  1.57  in  (40 

mm)  x  .25  in  (6.5  mm)  graphite/epoxy  plate  and  impactcr 
structure. 
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Figure  8  illustrates  the  time  history  of  axial  strain  which 
occurs  at  midspan  on  the  surface  of  the  plate  on  the  side 
immediately  opposite  the  impact  location.  Note  that  the  strain  is 
lower  when  transverse  shear  and  contact  deformations  are  included 
as  a  result  of  the  additional  energy  absorption  mechanisms  other 
than  flexure.  Table  2  presents  maximum  load  cell  force  and  maximum 
axial  strain  on  the  surface  opposite  the  impact  point  for  all 
plates  and  deformation  mechanisms  studied.  As  expected,  the  long 
plates  show  little  effect  of  neglecting  contact  and  transverse 
shear  deformations. 

As  previously  discussed,  finite  element  analysis  was  performed 
for  the  most  part  under  an  impactor  velocity  of  315  in/s  (8  m/s) 
and  no  gravity  effects.  Tests,  however,  were  carried  out  under 
velocities  which  varied  from  19  in/s  (480  mm/s)  to  170  in/s  (4300 
mm/s).  In  order  to  compare  finite  element  analysis  results  with 
test  results,  a  non-dimensionalization  of  analytical  data  from 
finite  element  runs  was  undertaken.  Using  a  one-degree-of -freedom 
system  model  of  a  mass  m  impacting  a  massless  spring  of  spring  rate 
k  with  initial  velocity  vq  in  a  gravity  field  (gravitational 
acceleration  g)  gives,  for  the  contact  force  P  between 
mass  and  spring 

P  =  vq (km) l/2sin(wt )  +  mg[l  -  cos(wt)]  (19a) 

where 

w  =  (k/m)1/2  (19b) 

and  t  is  time.  If  the  spring  is  a  long  plate  of  length  L  clamped 
at  either  end  with  thickness  h,  cross-sectional  moment  of  inertia 
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Table  2.  Maximum  impactor  load  cell  force,  P,  and  axial 

bottom-surface  strain,  e,  predicted  during  impact  of 
composite  plates  by  8.4  lb  (3.8  kg)  mass  at  315  in/ s  (8.0 


Plate 

ID 


Flexure  Flexure  Flexure  Flexure  Flexure  Flexure 
Contact  Sc  Shear  Only  Contact  &  Shear  Only 

Sc  Shear  &  Shear 


4.35 

4.23 

9.9 

8.16 

26.25 

20.2 

67.5 

56.0 


4.40 

8.25* 

10.0 

8.25 

28.7 

21.2 

70.0 

56.5 


NDG 

NDG 

10.1 

8.30 

31.5 

22.0 

72.0 

57. 


30.5 
60.0 
42.0 

60.5 
30.2 
46.0 
30.0 

49.5 


30.6 
81.0* 
42.5 
60.0 
37.0 

51.7 

33.7 
54.0 


NDG 

NDG 

43.0 

60.5 
38.0 

53.5 
37.0 
58.0 


NDG:  identifies  unacceptable  finite  element  results. 
*:  probably  incorrect  due  to  numerical  instability 


I,  and  material  plane  strain  extensional  stiffness  C^;  and  if  it 
is  further  assumed  that  the  static  flexural  deflection  curve 
represents  dynamic  deflection,  then 


k  =  192CXXI/L3 


(19c) 


The  first  term  in  equation  (19a)  is  the  force  due  to  impactor 
initial  velocity  and  the  second  term  is  the  effect  of  gravity.  A 
nondimensional  load  cell  force  p*FE  is  defined  for  finite  element 
results  which  adds  the  gravitational  force  as  follows: 


P*FE  =  PFE/tv0(km)1/2J  +  [g(km)1/2/(VQk) )[1  -  cos(wt)] 


(20) 


\\ 


v  '  V” 
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where  PEE  is  value  of  load  cell  force  from  finite  element  analysis 


ac  any  given  time  t.  The  flexural  strain  e  in  a  plate  undergoing 


dynamic  impact  loading  represented  by  equations  (19)  is  found  to  be 


e  =  ( 12h/L3 ) { ( vq/w) sin( wt )  +  (mg/k)[l  -  cos(wt)]} 


In  a  manner  similar  to  load  nondimensionalization,  a  nondimens ional 


strain  e*FE  which  adds  gravitational  effects  is  defined  as 


e*FE  =  eFE/[12v0h/(L2w) ]  +  [mgL2w/ ( 12v0hk) ] [ 1-  cos(wt)]  (22) 


where  eEE  is  the  axial  normal  strain  obtained  from  the  finite 


element  analysis.  Time  is  nondimensionalized  by  the  time  for  which 


a  mass  remains  in  contact  with  a  spring  for  a  simple  mass-spring 


system: 


t*  =  t/(7Tw) 


Nondimensional  load  cell  force,  p*FE,  and  nondimensional  strain 


e*FE  opposite  the  impact  location  on  the  plate  surface  have  been 


calculated  for  fully  elastic  finite  element  runs  of  all  plates  (II 


through  18).  A  complete  set  of  graphs  of  nondimensional  force  and 


strain  versus  nondimensional  time  are  presented  in  Appendix  B. 


C.  COMPARISON  OF  FLEXURAL,  CONTACT,  AND  SHEAR  DEFORMATION  EFFECTS 


Table  2  presents  load  cell  force  and  axial  normal  strain  on  the 


plate  surface  opposite  the  impact  point  obtained  from  finite 


element  analysis  with  (a)  all  elastic  deformations,  (b)  flexure  and 
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transverse  shear  deformations  but  no  contact  deformations,  and  (c) 
flexure  only  with  no  transverse  shear  or  contact  deformations.  As 
expected,  plates  with  large  aspect  ratios  ( length-to-thickness )  do 
not  appear  to  be  greatly  affected  by  neglecting  contact  or 
transverse  shear  deformations.  Results  for  plates  with  aspect 
ratios  less  than  40,  however,  show  that  both  contact  and  transverse 
shear  deformations  can  be  important. 

For  quantitative  comparison  purposes,  the  data  of  Table  2  has 
been  recast  into  Table  3  which  shows  the  percent  increase  in 
predicted  load  or  strain  which  occurs  when  ( 1 )  contact  deformations 
are  ignored  and  (2)  when  both  contact  and  transverse  shear 
deformations  are  ignored.  Differences  of  three  percent  or  less 
between  analysis  methods  were  not  considered  significant  and  are 
indicated  by  a  dash  in  the  table. 

Table  3  shows  some  results  which  were  unexpected.  It  was 
anticipated  that  the  shorter  plates  (smaller  aspect  ratio)  would 
show  the  greatest  differences  when  contact  and  shear  deformations 
were  neglected.  Instead,  the  shortest  plate  exhibits  less 
difference  from  the  fully  elastic  analysis  than  the  next  longer 
one.  Time  histories  of  the  finite  element  results  were  examined, 
and  it  was  found  that  the  major  reasons  for  the  larger  differences 
predicted  in  forces  and  strains  for  the  longer  15  plate  than  for 
the  shorter  17  plate  was  the  superposition  of  higher  order 
vibrational  harmonics  upon  the  first  mode.  It  appears  that 
artificially  stiffening  a  structure  by  neglecting  certain 
deformation  modes  can  create  natural  frequencies  which,  if  excited, 
may  override  the  fundamental  flexural  deformation  mode  giving 
larger  forces  and  strains  than  would  ordinarily  occur.  This  effect 
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Table  3.  Percent  increase  in  predicted  load  cell  force  and  axial 

normal  strain  caused  by  neglecting  contact  and  transverse 
shear  deformations  in  elastic  finite  element  analysis. 


Total 

%  Increase,  Load 
Ceil  Force 

Deformation  Mode 
Neglected 

Total 

%  Increase,  Axial 
Normal  Strain 

Deformation  Mode 
Neglected 

Plate 

Aspect 

Mtl . 

Contact 

Contact 

Contact 

Contact 

Ratio 

Angle 

(deg) 

Only 

&  Shear 

Only 

&  Shear 

11 

112 

0 

- 

NDG 

- 

NDG 

12 

112 

90 

+  95* 

NDG 

+  35* 

NDG 

13 

42 

0 

- 

- 

- 

- 

14 

42 

90 

- 

- 

- 

- 

15 

21 

0 

+9 

+  20 

+  23 

+  26 

16 

21 

90 

+  5 

+9 

+  12 

+  16 

17 

10 

0 

+4 

+  12 

+7 

+  23 

18 

10 

90 

- 

+9 

- 

+  17 

indicates  less  than  3%  difference  from  complete  elastic 

analysis . 


NDG:  identifies  unacceptable  finite  element  results. 
*  probably  incorrect  due  to  numerical  instability. 


is  not  expected  to  occur  for  every  plate  since  the  excitation  of 
higher  order  harmonics  will  depend  upon  plate  stiffness  (material. 


thickness,  and  length),  the  impactor  mass  and  stiffness,  and  the 
initial  impact  velocity.  Therefore,  one  set  of  impactor  and 
velocity  conditions  may  excite  one  structure  and  not  another.  This 
conclusion  is  important  because  it  means  that  "rules  of  thumb"  for 


deciding  when  to  and  when  not  to  use  a 
may  be  extremely  difficult  to  generate 
material  stiffness,  impactor  geometry, 
and  impactor  velocity  all  play  a  role. 


given  analysis  capability 
since  plate  geometry,  plate 
impactor  material  stiffness. 


§ 
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However,  it  is  clear  from  the  data  presented  in  Tables  2  and  3 
that  for  AS-3501  graphite  epoxy  plates  clamped  on  two  opposite 
ends,  elastic  contact  and  transverse  shear  deformations  play  an 
important  part  in  predicting  accurate  stresses  and  strains. 

Analyses  of  plates  with  length- to- thickness  ratios  greater  than 
40:1  do  not  require  the  inclusion  of  these  effects  to  obtain 
accurate  flexural  strain  predictions,  while  plates  with  20:1  and 
smaller  ratios  will  require  these  effects  if  less  than  3%  error  is 
desired.  It  is  apparent  that  each  effect  -  contact  and 
shear  -  accounts  for  anywhere  between  5%  and  20%,  and  they  both 
appear  to  be  equally  important  for  the  structures,  materials,  and 
impactor  analyzed.  It  may  be  that  these  effects  will  be  more 
important  for  plates  restrained  on  all  four  sides  and  for  spherical 
rather  than  cylindrical  impactor  tip  geometry. 

Further  conclusions  concerning  the  ability  of  elastic  analysis 
to  model  impact  of  composite  plates  will  be  drawn  from  comparison 
between  analysis  and  test  results  in  the  succeeding  sections  of 
this  report. 
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IV.  IMPACT  TESTS  OF  AS- 3 501  GRAPHITE/ EPOXY  PLATES 


EQUIPMENT  AND  TEST  SETUP 


1 .  Test  Specimens 

Tests  were  designed  to  duplicate  the  conditions  of  the 
analysis.  AS-3501  graphite/epoxy  laminated  plates  were  fabricated 
at  NAVAIRDEVCEN  using  two  stacking  sequences:  24-ply 
[ ( ±4 5 / 0 2 ) 2/±45/0/90 ]s  and  48-ply  [( ±45/02 ) 2/±45/0/90 ] 23 •  Specimens 
II  through  14  were  cut  from  the  24-ply  laminate  and  specimens  15 
through  18  were  cut  from  the  48-ply  laminate.  Odd  numbered 
specimens  had  the  high-stiffness  direction  along  the  plate’s 
longitudinal  axis,  while  even  numbered  specimens  had  the 
low-stiffness  direction  along  the  plate's  axis.  Figure  9  gives 
specimen  dimensions. 


2.  Test  Apparatus 


Plate  impact  samples  were  clamped  at  either  end  using  one  of 
two  clamping  apparati.  The  shorter  specimens  (13  through  18)  were 
c-clamped  between  3/8-in-  (9.5-mm-)  thick  steel  plates  to  a  steel 
frame  which  was  designed  especially  to  fit  the  NAVAIRDEVCEN  impact 
tower  bed.  The  longest  II  and  12  plates  were  clamped  in  a  modified 
rig  normally  used  for  plates  supported  on  all  four  edges  but  which 
was  modified  for  the  two-edged  support  situation.  The  primary 
difference  between  the  two  clamping  fixtures  was  the  length  of 
specimen  clamped  between  two  steel  plates.  Specimens  13  through  18 
had  clamp  plates  2  in  (51  mm)  long  on  either  end,  while  the  II  and 
12  specimens  were  clamped  1/2  in  (13  mm)  on  either  end.  Also, 
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SPECIMEN 

L 

in/mm 

C 

in/mm 

h 

in/mm 

STACKING  SEQUENCE 

11 

14/355 

0.5/13 

.128/3.25 

[(±45/02)2/±45/0/90]s 

13 

5.31/135 

2.0/51 

It 

ft 

15 

tf 

If 

.256/6.5 

If 

17 

2.56/65 

If 

If 

If 

12 

14/355 

0.5/13 

.128/3.25 

[ (±45/02)2/±45/0/90]2s 

14 

5.31/135 

2.0/51 

tf 

ff 

16 

If 

II 

.256/6.5 

If 

18 

2.56/65 

II 

ff 

1! 

Figure  9.  Laminated  graphite/epoxy  impact  test  specimen  geometry. 
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considerably  greater  clamping  force  could  be  exerted  by  the 
c- clamps  on  the  shorter  specimens  than  could  be  obtained  with  two 
eccentric  1/8  in  bolts  used  by  the  rig  for  clamping  specimens  II 
and  12.  As  described  below,  the  lesser  clamping  force  on  the  II 
and  12  specimens  allowed  membrane  forces  to  pull  the  specimens  from 
the  supports  invalidating  the  tests. 

Plates  were  impacted  using  the  NAVAIRDEVCEN  instrumented  impact 
tower  with  automated  data  storage  and  reduction  capabilities.  The 
tower,  shown  in  Figure  10,  consists  of  two  cylindrical  impactor 
guides,  the  impactor  which  contains  a  load  cell  for  measuring 
impact  force,  a  mechanism  for  drop  height  control/ impactor  lift  and 
automatic  release,  a  test  bed  with  specimen  supports,  and 
instrumentation  for  measuring  velocity  of  impactor  at  the  point  of 
impact.  A  Nicolet  two-channel  digital  storage  scope  is  used  to 
store  more  than  4,000  data  points  per  timed  event. 

The  main  body  of  the  impactor  (see  upper  left  inset,  Figure  7) 
is  fabricated  from  steel  plate.  It  contains  provisions  for  adding 
weights  to  control  impactor  mass.  For  the  tests  conducted  here,  no 
added  weights  were  used,  and  the  total  impactor  mass  was  either  8.4 
lbm  (3.8  kg)  with  a  2,000  lb  (9  kN)  load  cell  or  9.8  lbm  (4.4  kg) 
with  a  10,000  lb  (44  kN)  load  cell.  The  impactor  is  designed  to 
use  interchangeable  tips.  For  the  tests  conducted  here,  two 
special  cylindrical  impactor  tips  were  designed  and  manufactured  - 
one  with  a  2  in  (51  mm)  diameter  impact  surface  and  one  with  a  0.25 
in  (6.4  mm)  diameter  impact  surface.  Both  were  1.625  in  (41.3  mm) 
wide  -  slightly  wider  than  the  specimens.  Nearly  all  of  the  tests 
were  conducted  with  the  2,000  lb  (9  kN)  load  cell  and  the  0.25  in 
(6.4  mm)  impactor  tip. 
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Data  were  analyzed  using  impact  computer  software  and  plotting 
routines  developed  especially  for  the  NAVAIRDEVCEN  instrumented 
impact  tower  facility  by  Mr.  L.  W.  Gause.  This  software  provides 
time  plots  of  impact  force,  impactor  displacement ,  specimen  strain, 
and  absorbed  energy  for  the  impact  event,  and  plots  of  impact  force 
and  absorbed  energy  versus  impactor  displacement.  It  also* 
calculates  velocity  at  impact  and  maximum  values  of  force, 
displacement,  strain,  and  other  quantities  of  interest. 

B .  TEST  PROGRAM 


Nondestructive  impact  tests  of  all  eight  configurations  of 
AS-3501  graphite/epoxy  plates  were  performed.  Each  test  specimen 
was  strain  gaged  so  that  time  histories  of  axial  strain  at  midspan 
on  the  back-surface  (surface  opposite  the  impact  site)  could  be 
measured  and  compared  with  analysis.  Both  back-surface  strain  and 
impactor  load  cell  force  were  recorded  during  the  impact  event  on 
the  Nicolet  digital  storage  scope. 

Drop  heights  varied  between  1.0  in  (25  mm)  and  5.0  in  (125  mm) 
resulting  in  impact  velocities  from  about  2  ft/s  (0.6  m/s)  to  5 
ft/s  (1.5  m/s).  Care  was  taken  to  maintain  back-surface  strain 
less  than  0.005  to  avoid  structural  cracks. 

One  specimen  of  each  size  and  material  orientation  was 
subjected  to  three  consecutive  tests  at  the  same  drop  height  to 
demonstrate  repeatability.  Time  histories  of  impact  events  were 
found  to  be  nearly  identical  ( impact  loads  and  structural  strains 
were  virtually  the  same  from  test-to-test ,  varying  less  than  1%  in 
magnitude  and  duplicating  higher  order  harmonics  as  small  as  0.02 
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of  the  total  impact  period) .  Due  to  the  excellent  repeatibility  of 
test  data,  two  impact  tests  were  deemed  sufficient  for  each  test 
condition:  the  first  to  measure  force  and  strain,  and  the  second  to 
measure  impact  velocity  and  to  verify  that  the  specimen  had  not 
been  damaged  (identical  traces  of  impact  force  with  time  were 
considered  to  show  that  insignificant  damage  had  occurred) . 

Appendix  C  lists  the  impact  tests  performed  on  specimen 
configurations  II  through  18. 
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C .  RESULTS 


Contact  load  and  back  surface  strain  histories  for  each  test 
were  stored  on  floppy  diskettes  by  the  Nicolet  storage  scope.  The 
computerized  data  reduction  system  was  utilized  to  obtain 
experimental  verification  of  impact  velocities;  to  calculate 
impactor  head  displacement  and  absorbed  energy;  and  to  generate 
time  plots  of  impactor  force,  plate  axial  strain,  displacement,  and 
absorbed  energy. 

Figures  11  through  13  show  the  reduced  data  for  the  17  specimen 
impacted  at  3.6  ft/s  (1.1  m/s).  Similar  plots  have  been  obtained 
for  all  specimens  and  a  complete  set  is  provided  in  Appendix  D. 

It  is  noted  that  the  longest  plates  (specimens  II  and  12  with 
aspect  ratios  100:1)  pulled  loose  from  the  supports  due  to  large 
deflections  creating  significant  membrane  force  effects,  and  the 
data  may  be  invalid  for  comparison  with  analytical  predictions. 

For  comparison  with  the  analytical  predictions,  load  cell 
force,  plate  axial  normal  strain,  and  contact  impact  duration  were 
nondimensionalized  in  the  manner  used  for  the  analytically 
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Figure  12.  Impactor  load  cell  force  and  absorbed  energy  versus  time  for 
specimen  17,  impactor  velocity  3.6  ft/s  (1.1  m/s). 
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Figure  13.  Impactor  load  cell  force  and  absorbed  energy  versus  impactor 
displacement  for  specimen  17,  impactor  velocity  3.6  ft/s  (1. 
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determined  quantities,  but  without  the  gravity  term  which  is  not 
needed  for  the  test  data: 

nondimensional  load  cell  force  from  test: 

P*EXP  =  ^EXP^  t  v0  ( ,<m) 1/72  ^ 

nondimensional  strain  opposite  impact  location  from  test: 

e*EXP  =  eEXp/ C 12vgh/ ( L^w) ] 

where  Pexp  anc*  eEXP  are  test  impactor  force  and  plate  strain, 
respectively,  k  is  plate  stiffness  given  by  equation  (19b),  v0  is 
impact  velocity,  m  is  impactor  mass,  h  is  plate  thickness,  and  L  is 
plate  length  between  supports.  [Nondimensional  time  t*  is  given  by 
equation  (23).] 

Nondimensional  time  histories  of  impactor  force  and  plate  axial 
normal  strain  are  presented  in  Figures  14  and  15  for  the  17 
specimen  impacted  at  3.6  ft/s  (1.1  m/s).  Complete  nondimensional 
results  are  presented  in  Appendix  E. 

Nondimensional  test  results  are  summarized  with  finite  element 
results  in  Table  4.  The  first  line  of  p*  and  e*  data  for  each 
specimen  are  amplitudes  of  the  first  mode  force  and  strain 
responses,  respectively.  The  second  line  gives  the  maximum  higher 
mode  amplitudes.  The  maximum  amplitude  of  force  or  strain  may  be 
found  by  adding  first  mode  and  higher  mode  amplitudes.  The 
nondimensional  contact  impact  duration  t*MAX  is  the  time  for  which 
the  impactor  remains  in  contact  with  the  plate  during  the  initial 
impact  event  nondimensionalized  by  equation  (23).  Data  for  the 
long  plate  II  and  12  specimens  are  included  for  completeness  even 
though  they  may  be  invalid  as  discussed  above. 
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Figure  15.  Nondimensional  axial  normal  plate  strain  versus  nondimensional 
time  for  specimen  17  ,  impactor  velocity  3.6  ft/s  (1.1  m/s). 
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Table  4.  Nondimensional  impactor  force  p*,  plate  axial  normal 

strain  e* ,  and  contact  impact  duration  obtained  by 

finite  element  analyses  (FE)  and  tests  (EXP). 


1ST  MODE  FORCE  AND  STRAINS 

IMPACT 

STRAIN 

+  HIGHER  MODE  RESPONSE 

DURATION 

RATIO# 

1 

2 

3 

4 

5 

6 

7 

8 

SPM 

NO. 

L/h 

P*EXP 

P*FE 

e*EXP 

e*FE 

fc*MAXEXP 

P*MAXFE 

m 

[ 

11 

112 

0.85 

1 

1.11 

n 

1.05 

0.85 

1.01  ]$ 

r 

+  .1 

mm 

mm 

] 

H 

I 

s 

T 

13 

42 

0.65 

0.66 

1.11 

0.89 

1.2 

0.94 

1.25 

G 

I 

H 

F 

F 

15 

21 

0.50 

0.58 

0.87 

0.83 

1.45 

1.0 

1.05 

M 

N 

+  .02 

+  .15 

A 

E 

T 

S 

17 

10 

0.38 

0.58 

0.75 

0.81 

2.19 

1.47 

0.93 

L 

S 

[ 

12 

112 

n 

1.4 

in 

0.95 

0.8 

1.17  ] 

i 

DEB 

mm 

+  .08 

EH 

] 

L 

0 

s 

T 

14 

42 

0.60 

0.68 

1.3 

0.88 

1.24 

1.48 

W 

I 

+  .07 

+  .14 

+  .06 

I 

F 

I  1 

M 

F 

16 

21 

0.53 

0.62 

0.99 

0.85 

1.43 

1.16 

A 

N 

+  .06 

+  .16 

T 

E 

L 

S 

18 

10 

0.39 

0.69 

0.75 

0.89 

1.86 

1.28 

0.84 

S 

@  First  mode:  amplitude  of  fundamental  frequency  force  and  strain 
response . 

Higher  mode:  amplitude  of  frequency  response  for  second  and 
higher  mode  frequencies  combined. 

#  Strain  ratio:  ratio  of  plate  axial  normal  strain  on  the  back 
surface  obtained  experimentally  to  that  predicted  by  finite 
element  analysis. 

$  Experimental  data  for  the  long  II  and  12  plates  are  suspect  due 
to  support  slippage  during  the  test. 

HIGH  and  LOW  MATL  STIFFNESS  indicate  that  odd  numbered  specimens 
have  approximately  [±45/02 ]§  layup  and  even  numbered  specimens  have 
approximately  [±45/902]s  tayup.  Refer  to  Figure  9. 


-50- 


-yvViV.v.v. 


NADC-871 06-60 


Recall  that  impact  test  results  are  nondimensionalized  with 
respect  to  maximum  values  predicted  by  an  elementary  one  degree  of 
freedom  strength  of  materials  analysis:  nondimensional  quantities 
of  1.0  indicate  that  experimental  results  and  one  degree  of  freedom 
results  coincide.  As  expected,  test  results  (columns  3,  5,  and  7) 
for  the  13  through  18  specimen  tests  are  generally  closer  to  the 
one  degree  of  freedom  predictions  for  longer,  thinner  plates.  It 
is  not  likely  that  tests  should  agree  with  the  one  degree  of 
freedom  predictions  even  for  long  plates,  since  the  analysis 
ignores  higher  degree  vibrational  modes. 

Contact  and  transverse  shear  deformations  effectively  add 
flexibility  to  a  structure,  and  these  effects  are  more  important 
for  shorter  plates.  Therefore,  the  decrease  in  nondimensional 
impactor  force  and  plate  axial  normal  strain  with  decrease  in 
aspect  ratio  is  expected,  as  is  the  increase  in  contact  impact  time 
with  decrease  in  aspect  ratio. 

Except  for  the  longest  II  and  12  specimens,  the  data  appear  to 
be  consistent  and  in  accord  with  fundamental  principles. 
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V.  COMPARISON  OF  IMPACT  ANALYSIS  WITH  TESTS 

Table  4  presents  nondimensional  impactor  force,  plate  axial 
strain,  and  contact  impact  time  for  analytical  elastic  finite 
element  results  (columns  4,  6,  and  8)  in  addition  to  experimentally 
measured  values  (columns  3,  5,  and  7).  Column  9  is  the  ratio  of 
the  experimentally  measured  plate  axial  normal  strain  to  that 
predicted  by  the  finite  element  analysis. 

It  is  apparent  from  Table  4  that  the  small  displacement  elastic 
analysis,  even  with  transverse  shear  and  contact  deformations,  may 
be  inadequate  to  predict  structural  strains  for  composite  plates. 
Experimental  strains  were  from  15%  lower  to  45%  higher  than 
predicted  by  the  finite  element  elastic  analysis. 

The  shortest  plates  (aspect  ratio  10:1)  exhibited  smaller  axial 
normal  strains  and  impactor  forces  than  predicted.  Transverse  shear 
and  contact  viscoelasticity  may  account  for  these  differences.  If 
so,  the  effects  may  be  more  pronounced  for  plates  supported  on  all 
edges  and  for  sharper  impact  indenter  radii  than  the  1/8-in  (3-mm) 
radius  used  for  both  test  and  analysis. 

Experimental  axial  strains  in  the  longer  plates  (aspect  ratios 
greater  than  20:1)  were  greater  than  predicted  by  analysis.  It  is 
doubtful  that  viscoelastic  effects  can  account  for  this  difference, 
since  they  would  absorb  more  energy  and  reduce  flexural  strains. 
Membrane  effects  would  stiffen  the  structure  and  reduce  flexural 
strains.  Although  additional  axial  strains  would  be  imparted  by  the 
membrane  forces,  it  is  not  anticipated  that  they  would  increase 
total  axial  strains  beyond  those  which  would  be  predicted  by  small 
displacement  flexural  analysis. 
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VI.  CONCLUSIONS 

The  major  conclusions  from  this  study  are  as  follows: 

1.  Flexural  viscoelasticity  is  not  an  important  energy  absorption 
mode  during  low-velocity  impact  of  graphite/epoxy  composite  plates. 
The  degree  of  flexural  viscoelastic  behavior  required  to  fit 
experimental  impact  data  of  graphite/epoxy  plates  is  several  orders 
of  magnitude  larger  than  that  available  in  the  material. 

2.  Finite  element  analysis  indicates  that  elastic  contact  and 
transverse  shear  deformations  are  important  to  accurate  prediction 
of  structural  stresses  and  strains  for  plates  clamped  on  opposing 
edges  with  length-to-thickness  ratios  of  20  and  less.  Results  of 
purely  flexural  analysis  of  plates  with  aspect  ratios  of  40  and 
greater  appears  not  to  differ  from  results  of  analyses  containing 
elastic  contact  and  transverse  shear  deformation  effects. 

3.  The  neglect  of  contact  and  transverse  shear  effects  can  alter 
the  structural  response  of  an  impacted  plate  by  creating  a  higher 
order  resonant  response  (or  eliminating  one)  and  may  produce  large 
errors  in  stress  and  strain  predictions.  It  is  therefore  difficult 
to  develop  quantitative  "rules  of  thumb"  for  errors  resulting  from 
analytical  approximations,  as  the  dynamic  response  will  depend  upon 
plate  geometry,  plate  material,  support  conditions,  impactor  mass, 
and  impactor  structural  geometry  and  stiffness.  However,  present 
analyses  indicate  that  errors  between  5%  and  20%  can  result  for 
short  plates  by  neglecting  either  contact  or  shear  deformations. 
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4.  Membrane  effects  are  definitely  important  for  plates  clamped  on 
two  opposing  edges  having  aspect  ratios  of  100  or  greater. 

5.  Comparison  between  analysis  and  test  shows  that  elastic  contact 
and  shear  deformation  modeling  still  overpredicts  strains  in 
impacted  graphite/epoxy  plates  with  aspect  ratios  of  10  or  less. 
Transverse  normal  and  shear  viscoelastic  behavior  might  explain  the 
difference.  It  is  recommended  that  an  analysis  capability  which 
includes  these  viscoelastic  effects  be  developed  and  used  to  study 
stresses  and  strains  in  plates  supported  on  all  sides.  It  is 
probable  that  fully  supported  plates  impacted  with  a  spherical 
indenter  will  exhibit  greater  discrepancies  than  the  plates  used  in 
this  study  which  were  supported  on  two  opposing  ends  and  impacted 
with  a  cylindrical  indenter. 

6.  Experimental  strains  in  long  plates  (aspect  ratios  of  20  and 
greater)  can  be  as  much  as  48%  higher  than  analytical  elastic 
predictions.  This  might  be  due  to  two  effects: 

(a)  membrane  forces  during  large  deflections,  and/or 

(b)  significant  structural  natural  frequency  alteration  due  to 
improper  material  modeling. 

This  points  to  the  need  for  inclusion  of  membrane  force  effects  (e. 
g. ,  references  [12,  13]  or  [17])  in  the  analysis  of  plates, 
reinforces  the  desirability  of  modeling  viscoelastic  contact  and 
shear  behavior  of  the  material,  and  shows  the  advisability  of 
conducting  a  study  which  will  evaluate  these  effects  for  plates 
supported  on  all  edges. 


-54- 


J 

V 


NADC-d  7  1  Ot>-oG 


7.  It  is  well-known  that  the  prediction  of  accurate  deflections 
does  not  insure  that  stresses  and  strains  are  accurately  predicted: 
since  displacements  are  integrals  of  strains  and  integration  is  a 
"smoothing"  process,  there  are  many  examples  of  analyses  which  may 
have  only  3%  error  in  displacement  but  as  much  as  20%  error  in 
strains  and  stresses.  In  a  similar  fashion,  the  existence  of  the 
well-known  St.  Venant  effect,  where  self -equilibrating 

stress  states  decay  rapidly  with  distance  from  the  disturbance, 
generally  indicates  that  high  locally-generated  stresses  and 
strains  will  reduce  rapidly  with  distance  from  the  cause. 

Therefore,  the  errors  in  axial  normal  strain  reported  here  are  most 
likely  indicative  of  much  larger  errors  in  stresses  and  strains  in 
the  immediate  vicinity  of  the  contact  impact  location. 

8 .  Progress  has  been  made  to  determine  which  analytical 
capabilities  are  necessary  to  accurately  analyze  the  impact 
response  of  composite  plates  for  design  purposes.  However,  the 
program  has  raised  new  questions  as  well  as  answered  old  ones. 
Continued  investigation  into  this  problem  is  necessary  until  the 
answers  are  obtained.  Only  then  can  the  composites  community  be 
satisfied  that  it  can  accurately  and  efficiently  analyze  not  only 
structural  stresses  and  strains  (those  not  in  the  impact  region  and 
not  directly  affected  by  contact  stresses),  but  also  stresses  in 
the  immediate  vicinity  of  the  contact  region. 
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appendix  a 


CREEP  TEST  DATA  FOR  [  45]  AS-3501  GRAPHITE/ EPOXY 
LAMINATES 


TlMfc,  sec 


2IN  +45Z-45  CREEP  125#  C  (MWflt) 

ilKA/W.  tui  -  lysM/'t 


45  CREEP  250*  C  (wflww A«-> 
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/flit  i*  v 


1/2  In  STATIC  RAMP  LOAD  1000SEC  200* 


STATIC  RAMP  LOAD  ,  4SLC,  200*^Ai  *-<*» 


S  rc 


5/9/85  T# 1 :  1"  0-45]  320  LB  CREEP,  0  MIN  *  CL  ct&fP 

S Ph,  a.  tr*c  q 
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l&OLofltD 


5/9/85  T#  1  s  1“  C  +-45  3  320  LB  CREEP,  10  MIN.  Met  cae* pi 


S'  S  2  W  O! 
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(A/qj  OPT.) 
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\8o  L.OA& 
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appendix  b 


nondimensional  impactor  force  p*  an?T3^i^mtte  element 

BACK-SURFACE  STRAIN  e*  VS.  TIME  t*  fr°MtJ^cE  ELEMENT 
ANALYSIS  OF  IMPACTED  AS-3501  COMPOSITE  PLATES 
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APPENDIX  C.  LIST  OF  PLATE  IMPACT  TESTS,  TEST  CONDITIONS,  AND 
MAGNETIC  DISK  STORAGE  DATA 
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TEST 

DATA 

LOAD 

STRAIN 

DROP 

DATA 

TRACK 

COMMENTS 

NO 

TAKEN 

CALIB 

CALIB 

HEIGHT 

DISK 

NO 

*  ★ 

( kN/div) 

(V/ue) 

( in, mm) 

MBI* 

171 

P,e 

0.8 

6.0 

2.0, 

51 

2 

1 

172 

P,e 

0.8 

6.0 

2.0, 

51 

2 

2 

Rerun 

173 

P,e 

0.8 

6.0 

2.5, 

76 

2 

3 

174 

P,e 

0.8 

6.0 

2.5, 

76 

2 

4 

Rerun 

175 

P,  V 

0.8 

6.0 

2.5, 

76 

2 

5 

176 

P,v 

0.8 

6.0 

5.0,127 

2 

6 

177 

P,e 

0.8 

6.0 

5.0,127 

2 

7 

178 

P,e 

0.8 

6.0 

2.0, 

51 

2 

8 

Rerun 

179 

P,v 

0.8 

6.0 

2.0, 

51 

3 

1 

Rerun 

1710 

P,v 

0.8 

6.0 

2.5, 

76 

3 

2 

Rerun 

1711 

P,e 

0.8 

6.0 

2.5, 

76 

3 

3 

Rerun 

1712 

P/e 

0.8 

6.0 

5.0,127 

3 

4 

Rerun 

1713 

P/V 

0.8 

6.0 

5.0,127 

3 

5 

Rerun 

181 

P,v 

0.8 

6.0 

2.0, 

51 

3 

6 

182 

P,e 

0.8 

6.0 

2.0, 

51 

3 

7 

183 

P,e 

0.8 

6.0 

4.0,102 

3 

8 

184 

P,  V 

0.8 

6.0 

4.0,102 

4 

1 

185 

P,  V 

0.8 

6.0 

4.0,102 

4 

2 

Rerun 

186 

P,  e 

0.8 

6.0 

4.0,102 

4 

3 

Rerun 

151 

P,e 

0.8 

6.0 

2.0, 

51 

4 

4 

152 

P,v 

0.8 

6.0 

2.0, 

51 

4 

5 

153 

P,v 

0.8 

6.0 

3.5, 

89 

4 

6 

154 

P,e 

0.8 

6.0 

3.5, 

89 

4 

7 

161 

P  /  e 

0.8 

6.0 

1.0, 

25 

4 

8 

162 

P,V 

0.2 

6.0 

1.0, 

25 

5 

1 

163 

P,v 

0.2 

6.0 

1.5, 

38 

5 

2 

164 

P,e 

0.2 

6.0 

1.5, 

38 

5 

3 

165 

P,e 

0.2 

6.0 

2.0, 

51 

5 

4 

166 

P,v 

0.2 

6.0 

2.0, 

51 

5 

5 

131 

P,v 

0.2 

6.0 

1.0, 

25 

5 

6 

132 

P,e 

0.2 

6.0 

1.0, 

25 

5 

7 

133 

P,e 

0.2 

6.0 

1.7, 

43 

5 

8 

134 

P/V 

0.2 

6.0 

1.7, 

43 

6 

1 

135 

P,v 

0.2 

6.0 

2.0, 

51 

6 

2 

136 

P,e 

0.2 

6.0 

2.0, 

51 

6 

3 

141 

P,e 

0.2 

6.0 

0.5, 

13 

6 

4 

142 

P/V 

0.2 

6.0 

0.5, 

13 

6 

5 

143 

P,v 

0.2 

6.0 

1.0, 

25 

6 

6 

144 

P/e 

0.2 

6.0 

1.0, 

25 

6 

7 

Ill 

P,e 

0.2 

6.0 

2.0, 

51 

6 

8 

*** 

112 

P/V 

0.2 

6.0 

2.0, 

51 

7 

1 

*** 

113 

P,v 

0.2 

6.0 

3.5, 

89 

7 

2 

*  *  * 

114 

P/e 

0.2 

6.0 

3.5, 

89 

7 

3 

*  *  * 

115 

P  /  e 

0,2 

6.0 

3.5, 

89 

7 

4 

Rerun*** 

116 

P/V 

0.2 

6.0 

3.5, 

89 

7 

5 

Rerun*** 

117 

P/V 

0.2 

6.0 

4.25, 

108 

7 

6 

*** 

118 

P/e 

0.2 

6.0 

4.25, 

108 

7 

7 

*** 

119 

P,e 

0.2 

6.0 

12.0, 

305 

7 

8 

Grp 

1110 

P,e 

0.2 

6.0 

12.0, 

305 

8 

1 

Rerun, Grp 

1111 

P/e 

0.2 

6.0 

36.0, 

914 

8 

2 

Fail, Grp 
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TEST 

DATA 

LOAD 

STRAIN 

DROP 

DATA 

TRACK 

COMMENTS 

NO 

TAKEN 

CALIB 

CALIB 

HEIGHT 

DISK 

NO 

*  * 

( kN/div) 

(V/ue) 

( in, mm) 

MBI* 

121 

P,e 

0.2 

6.0 

1.0,  25 

8 

3 

*  ★  ★ 

122 

P,e 

0.2 

6.0 

1.0,  25 

8 

4 

Rerun*** 

123 

P,v 

0.2 

6.0 

1.0,  25 

8 

5 

*  *  * 

124 

P,  V 

0.2 

6.0 

2.0,  51 

8 

6 

*  *  * 

125 

P ,  e 

0.2 

6.0 

2.0,  51 

8 

7 

*  *  * 

126 

P ,  e 

0.2 

6.0 

4.0,102 

8 

8 

*  *  * 

127 

P,v 

0.2 

6 .  o 

4.0,102 

9 

1 

Grp 

128 

P,e 

0.2 

6.0 

12.0,305 

9 

2 

Grp 

*  Data 

disks 

are  labelled 

"McLaughlin  Beam 

Impact 

(disk  no 

. )" 

**  p>e  _  impact  load  and  plate  strain.  P,v  -  impact  load  and 
impactor  initial  velocity. 

***  Data  for  these  specimens  suspect  due  to  support  slippage. 

Grp  -  Large  deflections  and  support  slippage  allowed  plate  to  be 
fully  dislodged  from  grips. 

Fail  -  Flexural  failure  observed  in  specimen. 
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LORD-  150  Lbf  MASS-  8.41  Lbro  HEIGHT-  1.00  Ft 

STRAIN-1  10235  Micro  Inch/lnch 

RADIUS-  .13  Inch  LOAD  SCALE-  . 2KN/DIV 
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APPENDIX  E.  NONDIMENS IONAL  PLOTS  OF  IMPACTOR  FORCE  p*  AND  PLATE 
.AXIAL  NORMAL  STRAIN  e*  VERSUS  TIME  t*  DURING 
GRAPHITE, EPOXY  PLATE  IMPACT  TESTS 
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INSTRUMENTED  IMPACT  TEST 


7.*  ■»)  b-  .181  Ft  <  48.0  mm)  Q-I4708E+5  Lb /Ft -2  (7.kjJBE+4  HH.) 


INSTRUMENTED  IMPACT  TEST 


INSTRUMENTED  IMPACT  TEST 

I82LNGC I6R) 


•+6.6313E-01  t»  M  .0940E+00 


INSTRUMENTED  IMPACT  TEST 

I8GLNGC I6B) 


INSTRUMENTED  IMPRCT  TEST 

154 ( I5G) 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPACT  TEST 

I32CI1H) 


000  - 


T 


Hh 

>V 


II 


INSTRUMENTED  IMPACT  TEST 


181  Ft  <  40.0  mm)  Q—16200E+3  Lb /Ft  ^2  <?.?S?E+4  HP.) 


O-16ZB0E+5  Lb/Ft^Z  C7.737E*4  MP*) 


INSTRUMENTED  IMPACT  TEST 

I4KI2F1) 
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INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPRCT  TEST 

144 (I2H) 


INSTRUMENTED  IMPACT  TEST 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPRCT  TEST 


INSTRUMENTED  IMPACT  TEST 
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DISTRIBUTION  LIST 

NON-GOVERNMENT  ACTIVITIES  (continued) 


NO  . 
COPI 


NORTHROP  AIRCRAFT  CORP.,  One  Northrop  Avenue,  Hawthorne,  CA  90250 
(Attn:  Dr.  M.  Ratwani ,  B.  Butler  and  R.  Whitehead). 

PURDUE  UNIVERSITY,  School  of  Aeronautics  and  Astronautics, 

West  Lafayette,  IN  47907 

(Attn:  Dr.  C.  T.  Sun).  ...... 

PROTOTYPE  DEVELOPMENT  ASSOCIATES,  INC.,  1560  Brookhollow  Drive 
Santa  Ana,  CA  92705 

(Attn:  E.  L.  Stanton).  ...... 

ROCKWELL  INTERNATIONAL,  Columbus,  OH  43216 

(Attn:  M.  Schweiger).  ...... 

ROCKWELL  INTERNATIONAL,  Los  Angeles,  CA  90009 

(Attn:  Dr.  Lackman).  ...... 

(Attn:  W.  O'Brien).  ...... 

ROCKWELL  INTERNATIONAL,  Tulsa,  OK  74151 

(Attn:  F.  Kaufman).  ...... 

ROHR  CORP.,  Riverside,  CA  92503 

(Attn:  Dr.  F.  Riel).  ...... 

SIKORSKY  AIRCRAFT,  Stratford,  CT  06622 

(Attn:  S.  Garbo).  ...... 

J.  P.  STEVENS  &  CO.,  INC.,  New  York,  NY  10036 

(Attn:  H.  I.  Shulock).  ...... 

TELEDYNE  RYAN  AERONAUTICAL  CO.,  San  Diego,  CA  92138 

(Attn:  R.  Long).  ....... 

UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE,  300  College  Park  Avenue, 
Dayton,  OH  45469 

(Attn:  Dr.  J.  Gallagher).  ..... 

UNIVERSITY  OF  DELAWARE,  Mechanics  &  Aerospace  Eng.  Dept., 

Evans  Hall,  Newark,  DE  19711 

(Attn:  Dr.  R.  B.  Pipes,  Dr.  J.  R.  Vinson  and  Dr.  D. Wilkins 
UNIVERSITY  OF  OKLAHOMA,  Norman,  OK  73019 

(Attn:  Dr.  C.  W.  Bert,  School  of  AMNE). 

UNIVERSITY  OF  WYOMING,  Laramie,  WY  82071 

(Attn:  Dr.  D.  Adams).  ...... 

VILLANOVA  UNIVERSITY,  Villanova,  PA  19085 

(Attn:  Dr.  P.  V.  McLaughlin).  .... 

VIRGINIA  POLYTECHNIC  INSTITUTE,  Blacksburg,  VA  24061 

(Attn:  Dr.  K.  Reifsnider).  ..... 

WASHINGTON  UNIVERSITY,  School  of  Engineering  and  Applied  Science, 
Materials  Research  Laboratory,  Campus  Box  1087, 

St.  Louis,  MO  63130 

(Attn :  T .  Hahn )..  .  .  .  .  .  . 
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DISTRIBUTION  LIST 

NON-GOVERNMENT  ACTIVITIES  (continued) 


GENERAL  ELECTRIC  CO.,  Philadelphia,  PA  19101 

(Attn:  A.  Garber,  C.  Zweben).  ..... 

GREAT  LAKES  CARBON  CORPORATION,  New  York.,  NY  10017 

(Attn:  W.  R.  Benn,  Manager,  Market  Development). 

GRUMMAN  CORPORATION,  South  Oyster  Bay  Rd.,  Bethpage,  NY  11714 

(Attn:  R.  Hadcock).  ....... 

(Attn:  S.  Dastin).  ....... 

HERCULES  AEROSPACE  DIVISION,  P.O.  Box  210,  Cumberland,  MD  21502 

(Attn:  Mr.  D.  Hug).  ....... 

HITCO ,  1600  West  135th  Street,  Gardena,  CA  90249 

(Attn:  N.  Myers).  ....... 

ITT  RESEARCH  INSTITUTE,  Chicago,  IL  60616 

(Attn:  K.  Hofar).  ....... 

KAMAN  AIRCRAFT  CORP . ,  Bloomfield,  CT  06002 

(Attn:  Technical  Library).  ...... 

LEHIGH  UNIVERSITY,  Bethlehem,  PA  18015 

(Attn:  Dr.  G.  C.  Sih).  ....... 

LEONARD  ASSOCIATES,  INC.,  6  East  Avenue,  Mt.  Carmel,  PA  17851 

(Attn:  Mr.  L.  Marchinski).  ...... 

LOCKHEED-CALIFORNIA  CO. ,  Burbank,  CA  91520 

(Attn:  E.  K.  Walker).  ....... 

(Attn:  A.  Vaughn).  ....... 

(Attn:  A.  James).  ....... 

LOCKREED-MISSILES  &  SPACE  CO.,  1111  Lockheed  Way,  Sunnyvale,  CA  94086 
(Attn:  J.  A.  Bailie).  ....... 

LOCKHEED-CALIFORNIA  CO.,  Rye  Canyon  Research  Laboratory, 

Burbank,  CA  91520 

(Attn:  D.  E.  Pettit).  ....... 

LOCKHEED-GEORGIA  CO. ,  Marietta,  GA  30063 

(Attn:  Technical  Information  Dept.,  72-34,  Zone  26).  . 

LTV  AEROSPACE  &  DEFENSE  CO.,  Vought  Missile  &  Advanced  Program 
Division,  P.O.  Box  225907,  Dallas,  TX  75265-0003 
(Attn:  R.  Knight).  ....... 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY,  Technology  Laboratory  for 

Advanced  Composite,  77  Massachuetts  Avenue,  Cambridge,  MA  02139 
(Attn:  Dr.  P.  A.  Lagace).  ...... 

MATERIALS  SCIENCES  CORP.,  Spring  House,  PA  19477 

(Attn:  Dr.  B.  W.  Rosen).  ...... 

McDONNELL-DOUGLAS  CORP.,  St.  Louis,  MO  63166 

(Attn:  K.  Stenberg,  R.  Garrett,  R.  Riley,  J.  Doerr).  . 
McDONNELL-DOUGLAS  CORP.  Long  Beach,  CA  90846 

(Attn:  J.  Palmer).  ....... 

MpDONNELL-DOUGLAS  HELICOPTER  CO.,  Culver  City,  CA  90230 

(Attn:  J.  K.  Sen,  Trailer  2002).  ..... 

McDONNELL-DOUGLAS  HELICOPTER  CO.,  5000E.  McDowell,  M/S  B337 
Mesa,  AZ  85205 

(Attn:  Steve  Guymon).  ....... 
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DISTRIBUTION  LIST 
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NO  . 

COP 

ANAMET  LABORATORIES,  LOO  Industrial  Hyw.  ,  San  Carlos,  CA  94070 

(Attn:  Dr-  R.  Arnold).  .......  1 

ALCOA  DEFENSE  SYSTEMS  CORP . ,  16761  Via  delCampo  Court, 

San  Diego,  CA  92127 

(Attn:  D.  Myers).  .......  1 

AVCO,  Specialty  Materials  Div. ,  2  Industrial  Avenue,  Lowell,  MA  01851 

(Attn:  Mr.  W.  F.  Grant).  ......  1 

BATTELLE  COLUMBUS  LABORATORIES,  Metals  and  Ceramics  Information  Center 

505  King  Avenue,  Columbus,  OH  43201.  .  .  .  .  1 

BEECH  AIRCRAFT  CORP.,  4130  Linden  Avenue,  Dayton,  OH  45432 

(Attn:  M.  B.  Goetz).  .......1 

BELL  AEROSPACE  COMPANY,  Buffalo,  NY  14240 

(Attn:  F.  M.  Anthony,  Zone  1-85).  .....  1 

BELL  HELICOPTER  CO.,  Fort  Worth,  TX  76101 

(Attn:  M.  K.  Stevenson).  ......  1 

BENDIX  PRODUCTS,  Aerospace  Division,  South  Bend,  IN  46619 

(Attn:  R.  V.  Cervelli).  ......  1 

BOEING  CO.,  P.  0.  Box  3707,  Seattle,  WA  98124 

(Attn:  J.  McCarty,  J.  Quinliven,  and  Dr.  R.  June).  .  .  3 

BOEING  CO.,  Vertol  Division,  P.0.  Box  16858,  Philadelphia,  PA  19143 

(Attn:  R.  L.  Pinckney).  ......  1 

(Attn:  D.  Hart).  ........  1 

(Attn:  C.  Albrecht).  .......  1 

BOEINGCO.,  Wichita,  KS  67277-7730.  .....  1 

(Attn:  J.  Avery).  .......  1 

(Attn:  R.  Waner).  .......  1 

CABOT  CORPORATION,  Billerica  Research  Center,  Billerica,  MA  01821.  1 

DEPARTMENT  OF  TRANSPORTATION,  Kendall  Square,  Cambridge,  MA  02142.  1 

(Attn:  Dr.  Ping  Tong,  DTS  76,  TSC).  .  .  .  .  1 

DREXEL  UNIVERSITY,  Philadelphia,  PA  19104 

(Attn:  Dr.  P.  C.  Chou).  ......'  1 

(Attn:  Dr.  A.  S.  D.  W ang )..  .  .  .  •  .  1 

E.  I.  DuPONT  COMPANY,  Textile  Fibers  Department,  Chestnut  Run  Location 
CR701,  Wilmington,  DE  19898 

(Attn:  V.  L.  Bertarelli).  ......  1 

FAIRCHILD  REPUBLIC  CO.,  Farmingdale,  L.I.,  NY  11735 

(Attn:  Mr.  Frank  Costa).  ......  1 

GEORGIA  INSTITUTE  OF  TECHNOLOGY,  Atlanta,  GA  30332 

(Attn:  (L.  Rehfield)).  .......  1 

GENERAL  DYNAMICS /CONVAIR,  San  Diego,  CA  92138 

(Attn :  Dr .  R.  Dunbar )..  .  .  .  .  .  .  1 

GENERAL  DYNAMICS,  Fort  Worth  Division,  PO  Box  748,  Fort  Worth,  TX  76101 
(Attn:  J.  A.  F an t)..  .  .  .  .  .  .  1 

(Attn:  Composite  Structures  Eng.  Dept.).  .  .  .  .  1 
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GOVERNMENT  ACTIVITIES  -  (continued) 

NO.  OF 
COPIES 


N A V SHI? RAN 0 C E N ,  Annapolis,  MD  21403 

(Attn:  H.  Edleste in,  Code  2870).  .  .  .  .  .  1 

NRL,  Washington,  D.C.  20375 

(Attn:  Dr.  I.  W'olock,  Code  6122;  Dr.  C.  I.  Chang, 

and  Dr.  R.  Badaliance).  .....  3 

N5WC,  WHITE  OAR  LABORATORY,  Silver  Spring,  MD  20910 

(Attn:  Dr.  J.  Goff,  Materials  Evaluation  Branch,  Code  R-34  .  1 

(Attn:  Dr.  J.  M.  Augl).  ......  2 

ONR,  800  N.  Quincy  St.,  Arlington,  VA  22217 

(Attn:  A.  Kushner  Code  432/A;  Y.  Rajapak.se,  Code  1132SM)  .  2 


0NT ,  800  N.  Quincy  Street,  Arlington,  VA  22217 

(Attn:  Cdr.  D.  Brown,  (OCNR-212).  .....  1 

PLASTEC,  Picatinny  Arsenal,  Dover,  NJ  07801 

(Attn:  H.  Pebly).  .......  1 

(Attn:  Librarian,  Code  DRDAR-SCM-0,  Bldg.  351-N).  .  .  1 

U.  S.  ARMY  MATERIALS  RESEARCH  LABS,  DRXMR-PL, 

Watertown,  MA  02171 

(Attn:  D.  Op  linger ) .  .  .  .  .  .  .  .  1 


U.  S.  ARMY  APPLIED  TECHNOLOGY  LABORATORY,  USARTL,  (AVRADC0M), 

Ft.  East  is,  VA  23604 

(Attn:  J.  Waller;  T.  Mazna).  .....  2 

U.  S.  ARMY  AIR  MOBILITY  R&D  LABORATORY,  Ft.  Eustis,  VA  23604 

(Attn.  H.  Reddick).  .......  1 

U.  S.  ARMY  R6T  LABORATORY  (AVRADCOM) ,  Ames  Research  Center, 

Moffet  Field,  CA  94033 

(Attn:  r.  Imraen,  DAVDL-AS-MS  207-5).  .  .  .  .  1 

U.  S.  NAVAL  ACADEMY,  Annapolis,  MD  21402 

(Attn:  Dr.  R.  D.  Jamison,  Mechanical  Engineering  Department).  1 

DAVID  TAYLOR  NAVAL  SHIP  RESEARCH  &  DEVELOPMENT  CENTER, 

Annapolis,  MD  21402 

(Attn:  E.  T.  Camponeschi,  Code  2844;  R.  Crane,  Code  2844).  .  2 

DAVID  TAYLOR  NAVAL  SHIP  R&D  CENTER 
Bethesda,  MD  20084 

(Attn:  A.  Macander,  Code  1720).  .....  1 


NAVAIRDEVCEN,  Warminster,  PA  18974 
(Attn:  Code  8131). 

(Attn:  Code  09L2)  . 
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AFNAL ,  WPAF3 , 

DH  -3-* 

33 

(Attn: 

FI3EC, 

Dr.  G.  S endec. cyj). 

. 

- 

(Attn: 

FIB/L. 

Kelly,  W.  Goesch,  C.  Ramsey). 

. 

3 

(Attn : 

FIBCA) . 

•  •  •  a 

» 

1 

(Attn: 

FIBE/Mr 

.  D.  Smith). 

. 

1 

(Attn: 

MLBM/Dr 

.  J.  Whitney,  M.  Knight). 

2 

(Attn: 

MLB/F. 

Cherry).  . 

. 

1 

(Attn: 

MBC/’ Reinhart ) .  .... 

.  . 

1 

(Attn: 

AFWAL/MLSE/S .  Fecheck). 

. 

1 

DEPARTMENT  OF 

THE  AIR 

FORCE,  Bldg.  410,  Bolling  Air 

Force  Base  , 

Washington,  D.C.  20332 

(Attn:  Dr.  M.  Salkind,  Dr.  Amos).  .....  2 

DEFENSE  TECHNICAL  INFORMATION  CENTER  (DTIC),  Bldg. #5,  Cameron  Station 


Alexandria,  VA  22314 

(Attn:  Administrator).  .......  2 

FAA ,  Washington,  D.C.  20591 

(Attn:  J.  R.  Soderquist,  AW-I03).  .....  1 

F.AA,  Technical  Center,  Atlantic  City,  NJ  08405 

L .  N'eri,  Code  ACT-330;  M.  Caiar'a,  Code  ACT -033).  .  2 

NASA,  Washington,  D.  C.  20546 

(Attn:  Airframes  Branch,  FS-120).  .....  1 

(Attn:  OAST/RM,  Dr.  D.  Mulville).  .....  1 

NASA,  George  C.  Marshall  Space  Flight  Center,  Huntsville,  AL  35812 

(Attn:  E.  E.  Engler,  S&E-ASTN-ES) .  .....  1 

(Attn:  R.  Schwinghamer ,  S&E-ASTN-M).  .  .  ,  .  1 

NASA,  Langley  Research  Center,  Hampton,  VA  23365 

(Attn:  Dr.  J.  R.  Davidson,  MS  L88E;  Dr.  J.  Starnes,  MS-190; 

Dr.  M.  Mikulus,  H.  Bohan,  and  Dr.  C.  P.  Blakenship 
MS  189M).  ........  5 

NASA,  Lewis  Research  Center,  Cleveland,  OH  44135 

(Attn:  Dr.  C.  Chamis ,  MS  49-6;  M.  Hershberg,  MS  49-6).  .  2 

NAVAIRSYSCOM,  Washington,  D.C.  20361 

(Attn:  AIR-00 D4) .  .......  1 

(Attn:  AIR-530).  ........  1 

(Attn:  AIR-5302D).  .......  1 

(Attn:  AIR-5302).  .......  1 

(Attn:  AIR-5302F ) .  .......  1 

(Attn:  AIR-53032D).  .......  1 

(Attn:  AIR-9318) .  .......  1 

NAVPGSCHL ,  Monterey,  CA  95940 

(Attn:  Prof.  R.  Ball,  Prof.  M.  H.  Bank,  Prof.  K.  Challenger).  3 

NAVSEASYSCOM,  Washington,  D.C.  20360 

(Attn:  C.  Zannis,  SEA-05R25)).  .....  1 

NAVSEC,  Arlington,  VA  20360 

(Attn:  NSEC-6 101E) .  .......  1 


